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We discuss a relativistic diffusion in the proper time in an approach of Schay and Dudley. We 
derive (Langevin) stochastic differential equations in various coordinates. We show that in some 
coordinates the stochastic differential equations become linear. We obtain momentum probability 
distribution in an explicit form. We discuss a relativistic particle diffusing in an external electro- 
magnetic field. We solve the Langevin equations in the case of parallel electric and magnetic fields. 
We derive a kinetic equation for the evolution of the probability distribution. We discuss drag terms 
leading to an equilibrium distribution. The relativistic analog of the Ornstein-Uhlenbeck process is 
not unique. We show that if the drag comes from a diffusion approximation to the master equation 
then its form is strongly restricted. The drag leading to the Tsallis equilibrium distribution satisfies 
this restriction whereas the one of the Juttner distribution does not. We show that any function of 
the relativistic energy can be the equilibrium distribution for a particle in a static electric field. A 
preliminary study of the time evolution with friction is presented. It is shown that the problem is 
equivalent to quantum mechanics of a particle moving on a hyperboloid with a potential determined 
by the drag. A relation to diffusions appearing in heavy ion collisions is briefly discussed. 

PACS numbers: 02.50.Ey,05.10.Gg,25.75.-q 
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I. INTRODUCTION 

There were many attempts to generalize the diffusion 
in a way respecting relativistic invariance and causality 
([1][2][3][4][5]; for a review and further references see [6] 
[7]). In this paper we develop the approach initiated by 
Schay [1] and Dudley [2]. The diffusion process which 
should be considered as a relativistic analog of the Brow- 
nian motion is uniqely defined by the requirement that 
this is the diffusion whose four momentum stays on the 
mass-shell. We discuss Ito stochastic differential equa- 
tion [8] defined as a perturbation of the relativistic dy- 
namics on the phase space. As an example we consider 
the motion in an electromagnetic field [9] (its general- 
ization, a motion in the Yang-Mills field [10] [11], could 
be treated in a similar way). Relativistic stochastic dy- 
namics preserving the particle's mass has no normaliz- 
able Lorentz invariant equilibrium measure. We discuss 
drags which lead to an equilibrium probability measure 
for a large time but violate the Lorentz invariance. A 
stochastic process with such a drag would be an analog 
of the Ornstein-Uhlenbeck process. Covariant drags de- 
scribing a relativistic particle in a medium moving with 
a velocity V are discussed in [4]. In such a case the 
Lorentz transformation of the friction is compensated by 
the transformation of the velocity V. 

A diffusion process can be considered as a relativistic 
approximation to more complex many particle processes. 
In particular, a motion of a heavy particle in an environ- 
ment of a gas of light particles, described in a Marko- 
vian approximation by the master equation [12], could 
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be approximated by a diffusion process [13]. Such an ap- 
proximation is applied in a description of the quark-gluon 
plasma [14][15][16][17][18][11], for an electron in a back- 
ground cosmic radiation [19] [20] or a particle moving in 
a fluctuating metric [21]. We show that an equilibrium 
distribution consistent with the diffusion approximation 
to the master equation is severely restricted. The Tsal- 
lis [22] distribution satisfies this restriction whereas the 
Juttner distribution [23] and quantum distributions do 
not. We discuss the form of the relativistic diffusion equa- 
tion and compare it to the relativistic diffusions discussed 
in heavy ion collisions. 

We explicitly work out solutions to the stochastic dy- 
namics and its transition function in various coordinate 
systems. We determine the momentum distribution of a 
relativistic particle in external electromagnetic fields. We 
discuss the time evolution with a friction leading to an 
equilibrium. We show that such a dynamics is equivalent 
to an imaginary time evolution of a quantum mechanical 
particle moving on the hyperboloid in a potential deter- 
mined by the drag. We expect that standard quantum 
mechanics methods can be applied for a detailed approx- 
imation of the diffusive evolution. The explicit formulas 
may be useful for a comparison of theoretical predictions 
with experimental results of ultra-relativistic collisions 
when a gas of relativistic particles is formed. 



II. RELATIVISTIC DYNAMICS 

We are interested in random perturbations of the 
dynamics of relativistic particles of mass m. On the 
Minkowski space the dynamics of a relativistic particle 
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is described by the equations [9] 
dx» 1 



the embedding of H in R 4 . We obtain (expressing dpo 
by dpk from eq.(3)) 



dr 



~dT 



-P' 



m 



= K» 



(1) 



(2) 



where /i = 0,1,2,3 and is a force. The four- 
momentum p{r) of a relativistic particle defines the mass 
by the relation 



9jk = S jk - p 2 p 3 p k 
Then, g = m 2 c 2 p^ 2 and 

g jk = pk + ( mc y2 p j p k 

Hence, 



(9) 
(10) 



p 



A H = d\ + dl + d 2 + {mcy 2 p>p k d 3 d k + (mc)- 2 3p k d k 
Po(t) 2 - P i(t) 2 - p 2 (r) 2 - P3 (t) 2 = m 2 c 2 (3) where k = 2 ' 3 and a j = W We prefer another choice 



Eq.(3) (together with the positivity of energy) says that 
the momenta stay on the upper half TL + (defined by 
Po > 0) of the four-dimensional hyperboloid H. If eq.(3) 
is satisfied then from cq.(l) it follows that r has the 
meaning of the proper time. If eq.(3) is to be true then 
the force K 11 must satisfy the subsidiary condition 



= 



(4) 



(we use the convention of a summation over repeated 
indices). We add a random force k^dr = dp^ to eq.(2) 
writing it in the form 



4V = K^dT + dp" 



(5) 



The diffusion p H \t) on the hyperboloid H+ is uniquely 
defined. It is generated by the Laplace-Beltrami operator 
on H 



Ah = ■l F d j j k y/gd k 
V9 



(6) 



Here, g = det (gjk) and g^ is the metric on TL. If 
we define the expectation value over the sample paths 
(starting from (x,p)) of the diffusion process <j) T (x,p) — 
E[<P(x(t),p(t))\ (we denote the expectation values by 
E [...]) then (f> T satisfies the diffusion equation 



V m ox^ op^ 2 



(7) 



with the initial condition (f>. -y 2 has the meaning of a 
diffusion constant. We write it in the form 



7 = men 

Then, kT 2 has the dimension of time. 



(8) 



III. 



COORDINATES ON H 



A proper choice of coordinates may be useful for a 
solution of differential equations. The momenta (3) p 
(with po = \/p 2 + m 2 c 2 ) could be used as coordinates 
on 7Y + . Then, the metric tensor can be obtained from 



of coordinates (p + ,p a ) (where a = 1,2) 

P+ = PO + P3 

In such a case g ++ = p+ 2 (p a p a + m 2 c 2 ), g +a = -p a p^ 
and g ab = 5 ab . Then, g = m 2 c 2 p+ 2 , g ++ = (mc)~ 2 p\, 
g+ a = (w,c)- 2 p+p a , g aa = 1 + (m,c)- 2 p 2 a , g 12 = 
(mc)~ 2 pip2- Hence, 

A H = 81 + 81 + (mc)- 2 p\d\ + (mc)- 2 p\d\ 

+ (mc)- 2 pi,dl + 2(mc)- 2 p + p a d a d+ + 2(mc)- 2 p 1 p 2 d 1 d 2 

+ (mc)- 2 'Sp + d + + (mcy 2 3p a da 

(12) 

where d a = The formula (12) can be rewritten in a 
Lorentz invariant form 

Ah = -3(mc)"V^ + (mc)- 2 ( P y - ^ v p 2 )d lx d v 

(13) 

where rf v is the Minkowski metric. In order to derive 
eq.(12) from eq.(13) assume that p 2 = m 2 c 2 and the func- 
tion <j> in cq.(7) is expressed as a function of p + ,pi,p 2 . 

It is instructive to compare p with some other widely 
used coordinates (a change of coordinates in the relativis- 
tic diffusion equation is also discussed in [6] [24]). First, 
let us consider the analogues of spherical coordinates 

Po = mccosha (14) 

pi — mc sinh a cos 4> sin 6 

p 2 = mc sinh a sin <j) sin 

Ps = mc sinh a cos 6 
In these coordinates the metric is expressed as 

ds 2 = (mcfida 2 + (sinha) 2 ^) (15) 

where 

ds\ = d6 2 + (sinO) 2 d<jy 2 (16) 

is the metric on the sphere S 2 . The Laplace-Beltrami 
operator reads 

(mc) 2 A H — (sinha)" 2 9 Q ,(sinha) 2 i9 a + (smlic^^A^ 

(17) 
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where As 2 is the Laplace-Bcltrami operator on the 
sphere. Next, let us consider the Poincare coordinates 
(<7ij 92, 93) on H+ which are related to momenta p as fol- 
lows 



and 



mc 

P3 + Pa = — 

93 



( 2 2 2 \ 

P3 - Po = (9i + 9 2 + 9 3 ) 

93 



(18) 



Pi = 



mcqi 
93 



P2 = 



mcq 2 
93 



where (73 > 0. Then, the metric is 

ds 2 = (mc) 2 q 3 ~ 2 (dq 2 + dq\ + dq 2 ) 

and 

(mc) 2 A H = q 2 (d 2 + d\ + d\) - q 3 d 3 
IV. STOCHASTIC EQUATIONS 



(19) 



(20) 



A solution of the diffusion equation as well as expec- 
tation values of observables can be expressed by an ex- 
pectation value over the solution of stochastic equations 
[8]. We have discussed stochastic equations correspond- 
ing to the diffusion on Tt + and their solutions in [25] [26] . 
The stochastic equations on H+ are also discussed in 
[6] [27] [24]. We solve these equations in the case of the 
free motion (K = 0) using the Poincare coordinates or 
the light-cone coordinates. In the Poincare coordinates 
the diffusion process is a solution of the linear stochastic 
differential equations 



dq a = nq 3 db a 



(21) 



a = 1,2. 



dq 3 = - —q 3 dr + nq 3 db 3 = -K 2 q 3 dr + nq 3 o db 3 

where Stratonovitch differentials are denoted by a circle 
and the Ito stochastic differentials without the circle (the 
notation is the same as in [8]). The Brownian motion 
appearing on the rhs of eqs.(21) is defined as the Gaussian 
process with the covariance 



E[b a (r)b c (s)} = 5 ac min(T, s) 
The solution of eqs.(21) is 

q 3 (r) = exp(-K 2 r + Kb 3 (r))q 3 



(22) 



q a {T) = q a + K [ q 3 (s)db a (s) (24) 
Jo 

for a = 1, 2. The solution could be applied for a calcula- 
tion of correlation functions and the transition function. 
The transition function P T of the diffusion is a solution 
of the equation 



d T P = yA ff P 



(25) 



with the initial condition P (q,q') = g~^S(q — q'). We 
have [8] (we calculated the transition function from the 
solution of the stochastic equations in [26]) 

P T (a) = (27r« 2 r)-la(sinh ( 7)- 1 exp(-^ - ^) 

(26) 

where the geodesic distance a in the Poincare coordinates 
can be expressed in the form 



cosh a = 1 + (2q 3 q / 3 )- 1 ((q 1 - q[) 2 + (q 2 - q' 2 ) 2 + (q 3 - q' 3 ) 2 ) 

(27) 

Using eqs.(18) we can derive the differentials dp^ . In the 
light-cone coordinates 



we have 



p± = po ± p 3 



3k 2 



(28) 



dp + = k p+dr + np + o db + = —^-p + dr + np + db + (29) 

(here we denoted b 3 by b + ) 
3 

dp a = 2 K Padr + np a db + + ^db a (30) 

where a = 1, 2 

p_ can be obtained from the formula 



Then 



Let 



p_ = (m 2 c 2 +p a p a )p- 1 (31) 



+ K{p _ _ 2jg£) db+ + ^_ pa Q dK 



MP) = mW))] (32) 



where p(r) is the solution of stochastic equations (29)- 
(30) with the initial condition p. Then, 



(33) 



(23) where is defined in eq.(12). 
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The spatial momenta which are useful for a physical 
interpretation of the relativistic diffusion lead to a non- 
linear (Ito) Langcvin equation 



where 



with 



dpi = -/sV<fr- + e£(p)d& n 



e n = + (Po - mc)(mc) x p Vp™ 



(34) 



(35) 



(36) 



These equations have been derived earlier in [6] . We can 
solve these equations by means of a change of coordinates 
(18) applying the solutions (23)-(24) or (29)-(30). 



V. PHASE SPACE EVOLUTION IN AN 
ELECTROMAGNETIC FIELD 

The evolution of coordinates can be obtained as an 
integral over the proper time 

1 r 

x^{t)=x i1 + — p^{s)ds (37) 
m Jo 

In an electromagnetic field the momentum satisfies the 
equation 



d PfJt = —F^p v dr + dp" 



(38) 



Here, by p H we denote the diffusion on the hyperboloid 
defined in eqs.(29)-(30). In general, we obtain non-linear 
stochastic differential equations from eq.(38) (because 
eq.(31) for p_ is non-linear in momentum). Eq.(38) is a 
linear equation if the equation for p + does not involve p- 
on the rhs. The only case which leads to linear stochastic 
differential equations describes constant parallel electric 
and magnetic fields (or a special case when one of them 
is zero). Let 



a = 



mc 



(39) 



Assume that the only components of F are F 12 = B and 
F 30 = E. In such a case eqs.(38) read 

3 

dpi — aBp 2 dr + -n 2 pxdr + npxdb + + ^dbi (40) 



dp 2 = —aBpidr + -n 2 p 2 dr + np 2 db + + ^db 2 (41) 



dp + = aEp + dr + K, 2 p + dr + np + o db + (42) 

It is clear that the linear equations (40)- (42) can explic- 
itly be solved. The solution of eq.(42) is an elementary 
function 



The environment of electromagnetic waves or (in a quan- 
tized form) photons can be another source of diffusion. 
Let in cq.(2) = (F^ + Q^ v )p v where Q is a Gaus- 
sian electromagnetic field (depending only on the proper 
time) with the covariance 



E[Q^ v (t)Q <tp (t')] = (VnaVvp - V^ p V^)S(t 



(44) 



Let p(r) be the solution of cq.(2) with an external 
(deterministic) electromagnetic field F and a random 
(or quantum) electromagnetic field Q. Then, 4> T {p) = 
E[(j)(p(T;p))] is the solution of eq.(7). 

In order to derive the non-relativistic limit we assume 
that 

p + = mc = const (45) 
in the stochastic equations (29)- (30). Then, 
d+ = d 3 

in the diffusion equation (33). In the non-relativistic limit 
eqs.(38) for a constant electromagnetic field become lin- 
ear. The solution is expressed by the Ornstein-Uhlenbeck 
process [28]. 

VI. THE MOMENTUM DISTRIBUTION 

We are interested in a distribution of momenta of par- 
ticles coming out from a gas formed after heavy ion colli- 
sions. For this purpose we express the transition function 
(26) by the momenta. The relativistic invariant formula 
reads 



cosh a = ^ (mc) 2 pp' 

In terms of the (p + ,pi,p 2 ) coordinates we have 

2 cosher = 2a = m~ 2 cr 2 p+ L p+ X 
((PiP+ - p'iP+f 
+ (P2P+ - P 2 P+) 2 
+m 2 c 2 p'f + m 2 c 2 p\ 



(46) 



(47) 



As a function of a the geodesic distance a has the form 

a = Ma + Va 2 - 1) (48) 
The time evolution in the momentum coordinates is 

MP) = T T (j>{p) = ( d»(p')P T (p,p')<j>(p') (49) 



where d\i = d 3 ppQ 1 mc is the relativistic invariant volume 
measure /i. We express eq.(49) in various coordinate sys- 
tems inserting the transition function (26) in eq.(49) with 
proper volume elements dfi. In the coordinates (14) the 
Ricmannian volume clement is 



P+(t) = p + cxp(aEr + k 2 t + k6 + (t)) 



(43) 



djjL = (mc) 3 dad6d(f> sinh 2 a sin 9 



(50) 
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In the Poincare coordinates 

rf/x = (mc) 3 dqidq2dq3q^ 3 (51) 
From the invariance under Lorentz transformations A 
cosha(p,p') — cosher(Ap, Ap'). (52) 
We can choose 

Ap= (po,0,0,p 3 ) (53) 
Let us define the rapidity y by 

pa ± p 3 = m T cexp(±y) (54) 

where 



2 2 22,2 22,2,2 

nirpC = m c + p T = m c +p 1 +p 2 



(55) 



The rapidity transforms in a simple way under the 
Lorentz boost (with the velocity v ) in the (0, 3) plane 



V = V + - 

c 

Then, in the frame where pr = 



a(p,p') 



y-y 



(56) 



(57) 



The rapidity is also closely related to the variable a in 
the coordinates (14). We have 

cosh a = cosh a cosh a' — sinh a sinh a' cos o 2 

where o 2 is the geodesic distance on the unit sphere. 
Hence, in the Lorentz frame (53) if 02 =0 then a = 
y — y' = a — a'. 

From eqs.(29) and (43) we can see that the process 
P+(t) is an exponential of a Gaussian process. Hence, its 
probability distribution should be the log-normal distri- 
bution. We could calculate it from the general formula 
(using the transition function (26)) 



<Pr(p+) = J dp[dp 2 dp' + p' + 1 P T (p,p')<p(p' + ) 
= fdp> + P^(p + ,p' + ) ( j>(p' + ) 



(58) 



However, it is easier to derive it directly from the solution 
(43). So, for the diffusion in the electric field (40)- (42) 
we obtain 



P^ + \p + ,p' + )E - (27tk 2 tY 



exp (_^( aS + K 2 ) 2__^ (ln( ^ )) 2 



^(p' + ) K ~ 2aE p+ 1 - K 



(59) 



VII. THE EQUILIBRIUM DISTRIBUTION 

Let us define the time evolution of an expectation value 
of an observable in a state p (a measure on the phase 
space) by 



dp T <i 



dp(f) T 



(60) 



We say that a measure v is the invariant measure for 
the diffusion process (see [8]) if the expectation value in 
eq.(60) is time- independent, i.e. 



/ 



dv(p,x)<f) T {p 1 x) = const 



(61) 



Assume that (in general) the diffusion equation reads 
d T (f>T = Q4>T 

where 



and 



dpi m dx^ 



(62) 



(63) 



is the generator of the deterministic flow in the coordi- 
nates (10) or 



(64) 



dp a dp+ m dx^ 
in the coordinates (12). Let us write 

dp T = dxd 3 p$ T (65) 

Then, from eq.(60) 



where in the coordinates (10) 



(66) 



and 



V9 



(68) 



Let us write the invariant measure in the form 

dv = d 3 pdxy/g$ R = d 3 pdx$ $R (69) 

Here, rf 3 p$o = d 3 p^/g is the (not normalizablc; g is cal- 
culated below eq.(9)) equilibrium measure for Ah, i.e. 



A H $o = 



(70) 



Then, the invariant measure v for the diffusion (62) is 
determined by the solution of the equation (obtained 
by differentiating eq.(61) over r) 



or 



g*$ $ R = 



g$ R = o 



(71) 



(72) 
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where 



we have 



g = ^A H -p (^-R j + ^)p^ (73) 

It can easily be seen that if the (non-zero) limit of p T (as 
t — ► oo) exists then 



lim dp T = dxd 3 pp 1 $_r 



(74) 



(in the weak sense of the convergence of measures). We 
can express eq.(72) as an evolution equation in time a; 

d $n, - ^mp^AnQR - m^Rj + Pj£j)po 

(75) 

If Ri docs not depend on x° then eqs.(66) and (72) 
may have the same time- independent solutions determin- 
ing the static equilibrium distribution <& B . In general, 
eq.(75) is the transport equation for When x° — > oo 
then <j>^ tends to the x°-independent equilibrium distri- 
bution $e 



lim = $i 



(76) 



solving both eqs.(66) and (75). 

Eq.(72) can be considered as an equation for the drag 
if the equilibrium measure &e is fixed. It is not possible 
to obtain the equilibrium measure which is normalizable, 
Lorentz invariant and at the same time concentrated on 
the mass-shell (p 2 = m 2 c 2 ). We give up the explicit 
Lorentz invariancc. If wc still require the rotation invari- 
ance then it is useful to work in the spherical coordinates 
(14). In these coordinates eq.(72) for &e reads (we re- 
strict ourselves to the momentum dependence of $_e) 



7 



&h$e =pod a (w<& E ) 



where 



Y = u{a) 



d_ 

da 



(77) 



(78) 



of eq.(62) in the coordinates (14) has only one component 
u. From eq.(77) we obtain 



u> = -7 2 <9 Q ln$ B 



(79) 



If- 



then 



$e = ex p(~ /3 c Po) = exp(— /3mc 2 cosh a) (80) 



uj = 7 mc 3 sinha 

2 ' 

For the Bosc-Einstcin distribution 



§E = ( CX P {j3mc 2 cosh a) — 1 



(81) 



(82) 



= — -7 2 mc 2 /3sinha!^l — cxp(— /3mc 2 cosha)^ (83) 



The drifts can be inserted (after a change of coordinates) 
into the diffusion equation (7) or the stochastic equations 
(29)-(30) in order to determine the diffusive dynamics. 

The equilibrium distribution (80) determines the dif- 
fusion generator in spherical coordinates 



7 „0„,-2 d „0„,2 d ^ t~ A 



VtP u irv u ir 



S 2 



±K 2 f3cp%£ 



(84) 

where u = |p| and p° = \/m 2 c 2 + u 2 . If the drift i?- 7 for 
$e (80) is derived from eq.(72) in the p coordinates (10) 
then we obtain 



5 = irAfl - ^i^dcpop 1 -^- 
2 2 op 3 



(85) 



VIII. DIFFUSION EQUATION AS AN 
APPROXIMATION TO MASTER EQUATION 

The diffusion equation (7) could be considered as an 
approximation to the dynamics of a heavy particle em- 
bedded in a gas of light particles. The kinetic equation 
describing the flow conservation under a Markovian scat- 
tering process reads [12] [13] (here t = ^) 

(d t +cp : >p 1 -^- f )p(p,x) 

= k 2 j d 3 k(w(p + k, k)p(p + k, x) - w(p, k)p(p, x)j 

(86) 

where / d 3 kw(p 7 k) = 1 and k 2 w(p,k) is the probabil- 
ity that in the unit time the momentum p of the heavy 
particle is changed to p — k through scattering on light 
particles (we could assume k — k but this is not neces- 
sary). The diffusion equation can be obtained by means 
of the Taylor expansion in k [13]. In such a case the dif- 
fusion coefficients can be calculated using the formulas 



d 3 kw(p,k)k : > = k 2 {k j ) 



^D ij = ^k 2 j d 3 kw<p k) /.■'/,-' 



(87) 



(88) 



Then, 



|£>« = Ik 2 J d 3 kw(p,k)(k l - (k i ))(k j - (k j )) 
+ ±k 2 {k t )(W) = + \k- 2 C l & 



It follows that if w(p, k) > then the matrix 

M ij = D ij _ (90) 

must be positive definite. 

The drift and diffusion coefficients have the simplest 
meaning in the coordinates (10). In these coordinates, 



7 



comparing eqs.(7), (11), (75) and (86) we obtain the dif- 
fusion equation for $^ = p with 



Po 2 



From eqs.(10) and (90) 



Po 



(92) 



must be a positive definite matrix ( the metric tensor g tJ 
is defined in eq.(10)). Let us assume that 



= exp(/(/3f> c)) 
In such a case from eq.(72) we obtain 



(93) 



(94) 



— 2 ^ mc 
Pa 



l\2 



Hence, 

Mjk = l 2 {5]k - PjPkP 

+PjPk"f 2 (mc)- 1 P Q 2 p- 2 (pl - mcK 2 £- 2 p 2 (§ + \poM' 

For Jiittner [23] as well as Bose-Einstein equilibrium dis- 
tributions the longitudinal term in eq.(95) becomes nega- 
tive at large momenta. Hence, at large energies the diffu- 
sion equation could not be a good approximation to the 
master equation. However, for the Tsallis distribution 
[22] 



<f> E (x) = (1 + (q-l)x) — 



(96) 



we have 



Ri = - 2^Po/?« 2 c(l + (q- lj/fcpo)" 1 (97) 

xf'(x) in eq.(94) is a bounded function. As a consequence 
the longitudinal term in cq.(95) is positive definite. In 
such a case we can apply the diffusion equation as an 
approximation of the master equation at high momenta 
as well. The Tsallis distribution appears in many mod- 
els ranging from turbulence [35] to heavy ion collisions 
[17] [36] [38]. 

We return to the proper time evolution of the relativis- 
tic diffusion in an electromagnetic field of sec. 5. We re- 
strict ourselves to the static case without a magnetic field. 
Then, the electromagnetic potential is A = (Aq, 0,0,0). 
We notice that Q in eq.(62) consists of two parts : the dif- 
fusion in momentum and the drag which is a sum of the 
dynamical piece and a dissipative one. It is easy to see 
that in the equation (72) for the equilibrium measure the 
dissipative and dynamical parts must separately vanish 
(this is a version of the fluctuation-dissipation theorem) . 
Hence, 

dv = dPxdPppo 1 exp(/(/3(cp + eA ))) (98) 



is the equilibrium measure in a time-independent electric 
field E = — V^4q- The drag term of eq.(62) reads 



1 



() 



(91) Y = -K 2 [3cp a pi f((3p c+peA ) ——di A — + 



d p^ d 



dpi 



'dp 



m dx^ 
(99) 

If we suggest that the master equation (86) is an equa- 
tion in the proper time r (instead of t = — ) then there 
would be no damping factor 22£ in the drift (91). In such 
a case the longitudinal term would become negative for 
large momenta. As a consequence the relativistic diffu- 
sion equation (7) could not be a consistent approximation 
at large energies to the proper time master equation . 



IX. EVOLUTION WITH FRICTION 

In the case of time-independent drags (as discussed 
in sec. 7) the search of equilibrium can equivalently be 
treated as a study of either proper time evolution of xq 
-independent densities $ or a^-independent solutions od 
the transport equation (75). Without friction the proper 
time evolution is expressed explicitly by eqs.(49) and 
(26). With friction we do not expect explicit solutions. 
We must rely on approximate methods or computer sim- 
ulations. There is no normalizable invariant measure for 
the relativistic diffusion generated by Ah because the 
process is fast growing as can be seen from the solution 
(43) (without the electric field) 



E[p\{t)] -p 2 ,cxp(4 K 2 r) 



(100) 



The friction is damping the time evolution. Without the 
noise term the proper time evolution is determined by 
the solution of the equation 



Tr = R(P) 



(101) 



generated by the flow Y (63). The evolution in x is 
defined by the drag of eq.(75) 



dp m 

dx po 



R( P ) 



(102) 



As an example, for the generator (85) (without the dif- 
fusion) we have 

p(x ) = cxp(— k 2 /3mcx )p 

The evolution in proper time (101) is also expressed by an 
elementary function which has the asymptotic behavior 
(coinciding with the one in the time — ) 



|p(r)| ~ |p(r)|cxp(-K 2 mc 2 /3r)) 



(103) 



The dynamical systems (101)-(102) have a trivial limit 
when time tends to infinity. The diffusive spread- 
ing makes it non-trivial. The evolution with friction 
is described by the semigroup exp(r(|A^ + Y)). A 
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rough approximation of this evolution as a product 
exp(riAfl-) exp(rY) can be expressed as a diffusion act- 
ing on the deterministic flow (101). Such an approxima- 
tion is reliable only for a small time. 

Before we propose an exact method to approach the 
time evolution with friction let us explain it in the well- 
known case of the Ornstein-Uhlenbeck (OU) process. 
The evolution is generated by 



Id 2 d 
= 2 ~d£ 



(104) 



Applying the invariant measure for the OU process we 
can transform the generator (104) into the Hamiltonian 
of the harmonic oscillator 

H osc = --_ +T ^__ = "«P(- ^ 2 )GoucM^ 2 ) 

(105) 

We apply the method to the diffusion with friction. We 
have for the generator (62) (with the friction (94)) 



-exp(/)eexp(-/) = -^-AH + V 



(106) 



where 



^=y(c/?)V(.f' 2 + .r) + ^>P0/' 



We have obtained a Hamiltonian for a particle moving 
on the hyperboloid (3) in a potential V. We could apply 
either Hamiltonian methods of quantum mechanics to 
this model [33] or functional integration. In the latter 
case, let <j) = cxp(— /)?/>, then we have the Feynman-Kac 
formula 

r (p,aO = exp(-/) J B[exp(-/ T y(pf )ds)V(p? ,x T )] 

(108) 

where is the stochastic process (34)-(36) on the hy- 
perboloid (3) starting at the point p. Clearly the solu- 
tion (108) can also be expressed in the form 

<Pt(p,x) = E[(f>(p T ,x T )} 

where p T is the solution of the equation 

dp T = Kdr + dpf (109) 

with the initial condition p. 

Let us note that if / is determined by the relativistic 
Boltzman equilibrium distribution (80) (or Bose- Einstein 
or Jiittner) then the potential V in eq.(107) is growing 
quadratically. Such a quadratic growth can substantially 
change the large time behavior of the solution of the dif- 
fusion equation. In the case of the Tsallis distribution the 
potential V is bounded. The expansion in V of eq.(108) 
(which coincides with the Dyson expansion of quantum 
mechanics) is convergent for arbitrarily large time. If 
the evolution of the momentum distribution tends to the 
Tsallis distribution then this should also be visible from 
the solution (108) even far from the equilibrium. 



X. DISCUSSION 

It is expected that the relativistic diffusion will appear 
in relativistic models of plasma. Some calculations based 
on the master equation (86) have been performed al- 
ready in the eighties [14] [15] (for a recent review see [29]). 
The scattering probabilities can be calculated in quan- 
tum field theory Another approach applies the Wigner 
function [30] for a description of the particle phase space 
evolution [11][18] [31]. 

The approach based on the master equation leads to 
the diffusion equation (75) for the probability density. 
We may write the diffusion part of eq.(75) (no friction) 
as 



= didjDa 



d t A t 



(110) 



where 



ij = 7 {Sij - PiPjp )—+^ — ( m ) 



and 



Ai =r 



,3 1 



2 p mc 



Pi 



(112) 



(107) The total drift (with the friction) is 



Q = n 2 Pl (I + lf3p cf'((3p c)) — (113) 

V A Z / po 

The authors [15] [17] [32] [29] write a general form of the 
diffusion equation without specifying the diffusion coef- 
ficients. It follows from our work that the diffusion co- 
efficients D l i are defined by the relativistic invariance. 
Then, the functions Bn and B± in eq.(ll) of [17] and 
(the analogs) B\ and Bq in eq.(25) of [29] are uniquely 
determined. Relativistic invariance determines also the 
coefficient A of the drift in cq.(10) of rcf.[17] as ^ if 
friction is switched off. The drag defining the friction is in 
one to one correspondence with the equilibrium measure 
as discussed already in [17]. 

We have shown that the equilibrium measures result- 
ing from the diffusion approximation to the master equa- 
tion (86) which are selected by the additivity of entropy 
(Jiittner or Bose- Einstein) do not lead to a diffusion equa- 
tion which could be valid at arbitrarily high energies. The 
stochastic process with the generator (85) could have 
been considered as a relativistic generalization of the 
Ornstein-Uhlenbeck (OU)(see [4] for another definition 
of the relativistic OU process). However, the relativistic 
counterpart of the Maxwell-Boltzman equilibrium mea- 
sure does not seem to be restricted to the Jiittner distri- 
bution ( sec the discussion in [37] [17] [34] [38]). For this 
reason all the stochastic processes with the drifts (113) 
could be considered as relativistic OU processes (they 
have the usual OU non-relativistic limit). 

The transport equation resulting from quantum field 
theory has an application to ultra-relativistic collisions 
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(in particular to heavy ion collisions). The time evo- 
lution of observables (denoted by <j>) or the probability 
density of the diffusion process (denoted by $) has been 
expressed in this paper by analytic formulas or by equa- 
tions which could be solved numerically. The results are 
determined by the relativistic invariance and the form 
of the invariant measure. In this way the experimen- 
tal results concerning the particle momentum distribu- 
tion (available from RHIC [39] and future experiments 
on LHC) can decide whether the particles coming out 
from relativistic collisions can be described as diffusing 
in a gas of light particles. A diffusion equation can also 
be treated as a tool for a phenomcnological description of 
the scattering data in heavy ion collisions (sec [32] [40]). 

Finally, let us mention that although there is only one 
diffusion on the hyperboloid (3), nevertheless, there are 
many Markov processes on this hyperboloid. We have a 



general Levy-Khintchin representation formula for pro- 
cesses with independent increments. Among these pro- 
cesses we distinguish here the fractional diffusion (with a 
stable probability distribution) defined by 

d T <t> = {-A H y<i> (in) 

where < 5 < 1. The soluble case 5 = \ with the 
transition function 

P T (a) = rV2n-i^(r 2 +a^K 2 (V^T^) 

(115) 

(where K v is the Bessel function of the third kind) shows 
characteristic features of the fractional diffusion. An ap- 
plication of the fractional diffusion in heavy ion collisions 
has been suggested recently in [41]. 
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